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y—{ Abstract 
> ■ 

As a consequence of gravitomagnetism, which is a fundamental weak-field prediction of general relativity 
and ubiquitous in gravitational phenomena, clocks show a difference in their proper periods when moving 

^ along identical orbits in opposite directions about a spinning mass. This time shift is induced by the rotation 
of the source and may be used to verify the existence of the terrestrial gravitomagnetic field by means of 

<0 orbiting clocks. A possible mission scenario is outlined with emphasis given to some of the major difficulties 
q which inevitably arise in connection with such a venture. 

U '. 

Introduction 
> 

^ Among the essential predictions of General Relativity, gravitomagnetism (and gravitational waves) still lack 
^ direct observational evidence. The only direct reference to date for the existence of gravitomagnetism is 
based on the orbital data analysis of Lageos I and II, which suggests that the orbital precession due to 
the spinning Earth corresponds with the Einsteinian prediction within ~ 20% (Ciufolini et ah, 1998). The 
first space experiment directly designed to detect the gravitomagnetic field of the Earth will be Gravity 
Probe B, which shall measure the Lense-Thirring precession of a couple of test gyros carried by a spacecraft 
along a near-Earth polar orbit. In the following we describe an alternative space experiment to observe the 
terrestrial gravitomagnetic field, based on the temporal structure induced by the rotating Earth, and discuss 
the influence of the static part of the terrestrial gravitational field on the outcome of such a measurement. 



THE CLOCK EFFECT 



In the weak field and low velocity limit, the field equations of General Relativity reduce to a set of Maxwell- 
type equations, thereby giving rise in this approximation to a number of phenomena similar to those known in 
electrodynamics, among them the well-known Lense-Thirring precession, which is the gravitational analogue 
of the precession of a spinning dipole in a magnetic field. However, the gravitomagnetic field not only makes 
a test gyro to precess about the field lines but also affects the motion of a test body resulting in a difference in 
the proper period for co- and counter-revolving particles. Again, this can be easily understood by comparison 
with electrodynamics: Consider a test charge moving along a circular orbit about a central charge. When a 
weak and homogeneous magnetic field with orientation perpendicular to the orbital plane is switched on, the 
test charge experiences an additional Lorentz force, either adding to or opposing the Coulomb attraction, 
depending on the direction of the motion of the particle and the magnetic field, respectively. The particle 



can, however, still circle along the same path as without the field, but now with a somewhat increased or 
decreased velocity; exactly the same is expected to happen in the gravitomagnetic case. The reason for this 
is the formal equivalence between the Lorentz force and the Coriolis force, i.e. Larmor's theorem and its 
gravitational analogue (Mashhoon, 1993). It follows that among two particles along identical but opposite 
orbits about a rotating mass, one will move faster than the other, resulting in a difference of their revolution 
times. Numerically, for circular equatorial orbits this time difference (called gravitomagnetic clock effect) 
to leading order is found to be (Cohen and Mashhoon, 1993, Gronwald et al., 1997) 

At ~ 4vr^2 ~ 2 x 10 ~ ?s ' W 

where J and M are the angular momentum and mass of the central body, respectively, and the numerical 
value results upon inserting the appropriate values for the Earth. In comparison, for Sun and Jupiter, 
which have the largest specific angular momenta among all solar system bodies, the clock effect amounts to 
Ars un ~ 7 x 1CT 5 s and Arj up i t er — 5 x 1CT 5 s, respectively. 

Within the PPN formalism, for semiconservative theories and ignoring the Whitehead term and any pre- 
ferred location effects, an additional factor of (1 + 7 + a\/4)/2 appears in Eq. (1) (Mashhoon et al, 2000), 
showing that the clock effect depends, as expected, on the same PPN parameter combination as the gravito- 
magnetic precession of a gyroscope (e.g. Will, 1995). 

It should be emphasized that Eq. (1) does not contain Newton's gravitational constant G; in fact, this is the 
reason for the large numerical value in (1) with regard to relativistic standards. Another remarkable char- 
acteristic of relation (1) is the independence of the clock effect from the distance of the orbiting satellites, a 
feature that is to some extent reminiscent of the Aharonov-Bohm effect. Finally, the positive sign in Eq. (1) 
indicates that the corotating satellite is slower than the counterrotating one, apparently in contradiction to 
the Machian idea of the dragging of inertial frames. More about these interesting peculiarities can be found 
in Mashhoon et al, 1999. 

GRAVITATIONAL PERTURBATIONS 

In principle, the validity of Eq. (1) can be tested by means of two clocks carried aboard two satellites 
along identical but opposite trajectories about the Earth, where each clock is read out when it has exactly 
covered an azimuthal angle of 2tt. For an orbit of 7000 km radius, the satellite travels ~ 0.04 mas in 200 ns, 
whence the accuracy in the determination of the azimuthal closure must be at least of the same order, which 
makes extremely high demands on the tracking facilities. However, assuming no systematic errors in the 
position measurements, these requirements become less stringent with the increasing number of revolutions 
because of the accumulation of the gravitomagnetic clock effect. Another source of errors is due to the 
gravitational and non-gravitational perturbations on the satellites (Lichtenegger et al, 2000). While the 
latter may be avoided by means of drag-free satellites, the former must be carefully modeled in order to 
extract the gravitomagnetic effect out of the clock data. In the following, we will concentrate only on the 
variation of the orbital period due to the static part of the Earth's gravitational field (an analysis of the 
dynamical part will be given elsewhere). 

Eq. (1) applies to the difference in the sidereal orbital period of two oppositely circling satellites. The 
angular velocity of a satellite with respect to a fixed system is the sum to + Ocosi + M (see Iorio, 2000), 
where Co and Q are the time rates of the argument of perigee and the line of nodes, respectively, Xi is the 
(perturbed) mean motion and i denotes the orbital inclination. Upon expansion of the terrestrial gravity 
field in terms of spherical harmonics, the perturbation equations for the angular orbital elements read (e.g. 
Seeber, 1989) 
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Here, n and m are degree and order of the expansion, while p and q are summation indices running from 
to n and — oo to oo, respectively. Further, a, e and n denote the semi-major axis, eccentricity and the 
unperturbed mean motion of the satellites, and M B and R E are the mass and mean radius of the Earth. 
Finally, S nrnpq (u>, M) is a periodic function of the angular coordinates and F nmp (i) and G npq {e) represent 
the inclination and eccentricity functions, respectively, where a prime indicates differentiation with respect 
to the argument (Kaula, 1966). Secular perturbations are induced by the even zonal harmonics and affect 
the variables to, Q and M. while leaving a, e and i unchanged, i.e. keeping the right hand side of (2) 
constant. Because G npq (e) represents an infinite series in e with its lowest power equal to \q\, for circular 
orbits G npq (e = 0) = 5o q , where Sij is the Kronecker delta. Thus, to first order, the orbital period modified 
due to secular perturbations becomes in the limit e — > by means of Eqs. (2) 
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In the following, in accordance with Eq. (1), we will consider orbits with e = i = only. The most relevant 
perturbation by far is due to C20, i-e. the lowest zonal harmonic, resulting in a difference with respect to 
the unperturbed period of AP 2 s q c ~ —15.7 s for an orbit of 7000 km radius. Because the time dependence of 
Eqs. (2) is contained in the function S, we can estimate an upper limit for the change in the revolution time 
due to periodic perturbations by replacing S with its amplitude {C\ m + £^ m ) 1//2 and again calculating the 
difference from the unperturbed motion. The first non-vanishing periodic contribution is due to n = m = 2, 
having an "amplitude" of |Ai-22 a:E | ~ 0.158 s with a period of the perturbation of -P^ er * — 52 min for a 
prograde orbit. In fact, only those combinations of n and m do not vanish for which the sum n + m is even. 
Based on NASA's EGM96 gravitational field model (Lemoine et ai, 1998), Table 1 lists all non- vanishing 
secular contributions to the orbital period as well as the amplitudes (in conformity with the above meaning) 
and periods of the sectorial and tesseral perturbations up to order and degree 6. 
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Table 1. Secular and periodic changes in the orbital period due to the gravitational field up to degree and 
order 6 for an orbital radius of 7000 km (the index + denotes a prograde orbit). 

According to Table 1, the perturbation in the orbital period of a satellite due to the nonsphericity of the 
Earth can be up to 8 orders of magnitude larger than the gravitomagnetic "perturbation". The accurate 
modeling of the gravitational perturbations of the Earth is limited by the uncertainty in the harmonic 
coefficients C nm and S nm . Figure la illustrates the mismodeling 6P of the period induced by the secular 
perturbations up to degree 100, based on the accuracy of the even zonal harmonic coefficients given in the 
EGM 96 model. The three lines correspond to an orbital radius of 7000, 8000 and 12000 km, respectively. 
As can be seen, for a near Earth orbit the accuracy of about the first hundred spherical harmonics is less 
than required, while for a 12000 km orbit only the accuracy of the coefficients up to degree ~ 20 are of 
importance. Figure lb shows the maximum mismodeling of the periodic perturbations due to the uncertainty 
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Figure 1: Mismodeling of the period for three different orbital heights: (a) error due to the even zonal 
perturbations C n o, (b) error due to the uncertainty in the coefficients C n 2 and S n 2- 

of the harmonic coefficients C n 2 and S n 2- For a 7000 km orbit, the error is up to 1000 times larger than the 
effect being studied and decreases to a factor of 100 at a distance of 12000 km. However, in spite of these 
large errors, the influence of the periodic perturbations is of less significance because they will cancel out 
as well as will be masked by the clock effect after a sufficient number of revolutions due to its accumulative 
character. 

SUMMARY 

The gravitomagnetic clock effect, which involves a coupling between the orbital motion of a satellite and 
the rotation of the Earth, can be considered as an interesting alternative way to verify the existence of the 
gravitomagnetic field predicted by General Relativity. We have examined the influence of the static part of 
the gravitational field of the Earth on the period of a satellite and have estimated, based on the EGM96 
gravity field model, the expected mismodeling with respect to the observation of the gravitomagnetic clock 
effect. While periodic perturbations, which involve a mismodeling of 100-1000 times larger than required, 
may cancel out, secular perturbations must be handled carefully. Assuming an experimental error of ~ 10%, 
we find the residual uncertainty exceeding the required accuracy by a factor of 100. Thus, a successful obser- 
vation of the clock effect has to await a distinct improvement in the determination of the gravitational field of 
the Earth. This may be expected within a couple of years based on the data of the upcoming GOCE mission. 
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